Quasiparticle Energies and Band Gaps of Graphene Nanoribbons 
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We present calculations of the quasiparticle energies and band gaps of graphene nanoribbons 
(GNRs) carried out using a first-principles many-electron Green's function approach within the 
GW approximation. Because of the quasi-one-dimensional nature of a GNR, electron-electron inter- 
action effects due to the enhanced screened Coulomb interaction and confinement geometry greatly 
influence the quasiparticle band gap. Compared with previous tight-binding and density functional 
theory studies, our calculated quasiparticle band gaps show significant self-energy corrections for 
both armchair and zigzag GNRs, in the range of 0.5-3.0 eV for ribbons of width 2.4-0.4 nm. The 
quasiparticle band gaps found here suggest that use of GNRs for electronic device components in 
ambient conditions may be viable. 
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Graphene, a single atomic layer of graphite, has been 
successfully produced in experiment [H, 0, 0] , which has 
resulted in intensive investigations on graphene-based 
structures because of fundamental physics interests and 
promising applications 0, 0, H B B E3|- When 
graphene is etched or patterned along one specific direc- 
tion, a novel quasi-one-dimcnsional structure, a strip of 
graphene of nanometers in width, can be obtained which 
is referred to as a graphene nanoribbon (GNR). The 
GNRs are predicted to exhibit various remarkable prop- 
erties and may be a potential elementary structure for 
future carbon-based nanoelectronics 
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particular, as a fundamental factor in determining trans- 
port and optical properties, the electronic band structure 
of GNRs has been the subject of great interest. 

Depending on specific GNRs, previous studies using 
tight-binding or massless Dirac fermion equation ap- 
proaches have predicted GNRs to be either metals or 
semiconductors [H, El E3, El El 0] • Whereas, den- 
sity functional theory (DFT) calculation showed that 
all zigzag-edged and armchair-edged GNRs have a finite 
band gap when relaxation of the structure or spin polar- 
ization is considered H 21 1. Recent experiments have 
reported finite band gaps in all the GNRs that have been 
tested [H, 111- However, it is well established [II] that 
the Kohn-Sham eigenvalues from DFT calculation are 
inappropriate to describe the band gaps of semiconduc- 
tors. The disagreement between the Kohn-Sham band 
gap and experimental data is worse for nanostructures 
because of the enhanced electron-electron interaction in 
those systems. On the other hand, first-principles calcu- 
lation based on many-body perturbation theory, such as 
the GW approximation 2J, |25| , has been shown to be re- 
liable for obtaining quasiparticle band gaps of nano-sized 
semiconductors [2(3, l27|, |28|, |29( . Motivated by the impor- 
tance but the lack of accurate knowledge about quasi- 
particle band gaps of the GNRs and by the successes of 





FIG. 1: (Color online) (a) A ball-stick model for an 11-AGNR 
which has 11 C-C dimer lines making up its width w. Hy- 
drogen atoms (white balls) are used to passivate the edge 
CT-dangling bonds, x, y and z are the Cartesian coordinates, 
(b) A ball-stick model for a 6-ZGNR which has 6 zigzag chains 
along the z direction. 



the GW approximation for nano-size semiconductors, we 
carry out a first-principles calculation using the GW ap- 
proximation to determine the quasiparticle energy spec- 
trum and the band gaps of the GNRs. 

We consider two common types of GNRs. Their struc- 
tures are shown in Fig. [TJ The left one, called armchair 
GNR (AGNR), has armchair-shaped edges; the right one, 
called zigzag GNR (ZGNR), has zigzag-shaped edges. 
The dangling er-bonds at the edges are passivated by hy- 
drogen atoms. The structures of the GNRs studied here 
are fully relaxed according to the forces and stress on 
the atoms. Following conventional notation, a GNR is 
specified by the number of dimer lines or zigzag chains 
along the ribbon forming the width, for the AGNR and 
ZGNR respectively, as explained in Fig. [T] For example, 
the structure of Fig. [T] (a) is referred as a 11-AGNR and 
the structure in Fig. Q] (b) as a 6-ZGNR. In addition, 
when referring to the width of a GNR here, we define the 
width without including the hydrogen atoms at the edge, 
as shown in Fig. [T] 

Following the approach of Hybertsen and Louie [25j |. 
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we first obtain the electronic ground state with DFT 
within the local (spin) density approximation [L(S)DA]. 
Then, the quasiparticle energies are calculated within the 
GW approximation to the electron self energy. Norm- 
conserving pseudopotentials [30| and the plane-wave ba- 
sis are used. In this calculation, k-grid is sampled uni- 
formly along the 1-D Brillouin zone. To assure that the 
quasiparticle energies are converged to within 0.1 eV, a 
1x1x32 k-point sampling is used for AGNRs and a 1x1x64 
k-point sampling for ZGNRs. Since the supercell method 
is used in this calculation to mimic isolated GNRs, we use 
a truncated Coulomb interaction to eliminate the image 
effect between adjacent supercells 31, 32, 33|. Consider- 
ing the geometry of the ribbons, we employ a rectangular- 
shape Coulomb truncation as 

Vc = \ 0(|a:| - x c ) 6(\y\ - y c ) 9(\z\ - z c ), (1) 
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y + z 2 is the distance between two 



where r = \/ x 2 + y 2 
electrons; x c , y c and z c are cutoff parameters. As dis- 
cussed in previous studies 3l| , the dimension of the unit 
cell has to be 2x c x 2y c x 2z c . Because of the single layer 
structure of GNRs, the truncation lengths x c and z c are 
fixed for all GNRs in our calculations. The unit cell vol- 
ume is linearly dependent on the width of a GNR, and 
the number of plane waves needed is also scaled linearly 
with the width of ribbon, which significantly reduces the 
cost of the computation. 

Another important aspect of the calculation is that we 
have to include the spin degree of freedom to account for 
the spin polarization in ZGNRs. It is shown that the 
static polarizability matrix is diagonal in spin space [3] ■ 
Combining this with the fact that the bare Coulomb in- 
teraction is independent of the spin degree of freedom, 
the spin-polarized GW calculation proves to be almost 
the same as the non-spin-polarized case with the excep- 
tion of replacing the static polarizability with the sum 
of its two diagonal spin components. The details of the 
spin-polarized G W calculation can be found in Ref [34[ • 

The LDA and quasiparticle band gaps of eleven arm- 
chair GNRs are shown in Fig. [2j As is found in the LDA, 
the quasiparticle band structure has a direct band gap at 
the zone center for all AGNRs studied. In addition, the 
band gaps of the three families of n- AGNRs, which are 
classified according to whether n=3p + l,3p + 2 or 3p (n 
is the number of dimer chains as explained in Fig 1, and 
p is an integer) , show qualitatively the same hierarchy as 
those obtained in LDA (E^ +1 > Ef > Ef+ 2 ^ 0). 

However, the GW self-energy corrections to the band 
gap, Eg, are significant for all the AGNRs. The correc- 
tions are from 0.5 to 3 eV for the AGNRs in Fig. [2] with 
width from 1.6 to 0.4 nm, which are much larger than 
those found for bulk graphite or diamond [25| . A weaker 
screening contributes to this enhanced self-energy correc- 
tion because the GNRs are isolated and surrounded by 
vacuum that does not screen the Coulomb interaction. 



FIG. 2: (Color online) Variation of band gaps with the width 
of AGNRs. The three families of AGNRs are represented by 
different symbols. The values of the same family of AGNRs 
are connected by solid lines as guides to the eyes. The open 
symbols are LDA band gaps while the solid symbols are the 
corresponding quasiparticle band gaps. Dashed arrows are 
used to indicate the self-energy correction for the smallest 
width ribbon of each of the three families of AGNRs studied. 



In addition, the confined geometry (one-dimensional na- 
ture) of the GNRs enhances the effect of electron-electron 
interaction, which further enlarges the self-energy cor- 
rection. This kind of enhanced self-energy correction is 
also found in other nanostructures such as nanotubes and 
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nanowircs 



The band gaps from both the LDA and GW calcu- 
lations show clearly size dependence in Fig. [2] because 
of quantum confinement. Under a hard- wall boundary 
condition used in previous works, an inverse relation, 
Eg cx l/(w + Wq), is widely applied to characterize the 
size dependence of the band gap in AGNRs, where w is 
the width defined in Fig. [2] and wq is a small constant 
(2.4 A). This size-dependence of the band gaps describes 
the tight-binding and LDA results well [21]. However, 
the boundary condition for the GNRs is not strictly a 
hard wall condition, and the electron distribution will 
leak out of the boundary more or less. Therefore the ef- 
fective width of GNRs should be larger than the physical 
width w. Considering this effect, we use the formula 



E„ — r, 

w + Wq + o 



(2) 



to fit the band gap values in Fig. [2J and the fitted results 
are given in Table I. For LDA results, the parameter S 
is close to zero or a little bit negative. For GW results, 
the parameter S found is around 1.5 to 2.9 A. Therefore 
the correction to the effective width for each edge is only 
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TABLE I: Fitted parameters for the LDA and GW band gaps 
of AGNRs according to formula (2). The value of fitted 5 
reflects an effective width correction. 





LDA 




GW 




family 


a(eV ■ A) 


5(A) 


a(eV ■ A) 


5(A) 


3p+l 


15.8 


0.1 


44.4 


1.8 


3p+2 


3.0 


-0.4 


14.6 


1.3 


3p 


7.6 


-1.7 


31.3 


2.9 
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around. 0.75 to 1.45 A, which may reflect the non-hard- 
wall nature of the confinement. In the limit of wider 
GNRs, 8 and wq can be ignored, the trend reduces to 
Eg oc 1/w. 

Fig. 11 (a) shows the LSDA bandstructure of 12-ZGNR. 
There are two notable characteristics in the electronic 
structure of ZGNRs: 1) the top of valence band and 
the bottom of conduction band are composed of mainly 
edge states; and 2) the spin interaction introduces a fi- 
nite band gap in the ZGNRs. As shown in Fig. [3] (b) 
and (d), the self-energy corrections to the LSDA energy 
gaps in ZGNRs are similar with those in AGNRs, and the 
corrections enlarge the band gap by 0.8 eV to 1.5 eV for 
the ribbons studied. The spin polarization changes the 
screening type of ZGNRs from that of a metal to that of 
a semiconductor. Therefore a significant self-energy cor- 
rection is resulted as in the case of the AGNRs. Based on 
the discussion of AGNRs, we also try to fit the width de- 
pendence of the quasiparticle band gaps in Fig.[3](b). We 
fit the results directly with a functional form of l/(w+5). 
The fitted S of LSDA is almost zero, and it is 16 A for the 
GW values, which is much larger than that in AGNRs. 
This is not unexpected, because it is the spin interaction 
between electrons close to the edge that induces the finite 
band gap in ZGNRs. Therefore we do not expect a sim- 
ple quantum-confinement effect, a X/w size-dependence 
of the band gap, in such narrow ZGNRs. 

Unlike the band gap (A ) located around three-fourth 
of the way to the Brillouin zone edge (Fig. [3] (a)), the 
energy gap at the zone boundary (A 1 ) is not sensitive to 
the width of ZGNRs as seen in Fig.[3](c). Previous tight- 
binding calculations [33] show that the profile of edge 
states decays to the center of ZGNRs with the factor of 
e~ ar , where a = -^ln|2cos^| (?f < kc< n,c is the 
lattice constant of ZGNRs along the z direction). As a 
result, the band edge states close to the zone boundary 
are highly confined at the edge of ZGNRs. Because of 
their dominant edge-state character, these states are not 
sensitive to the width of the ribbons, hence the gap A 1 
is virtually independent of width. 

Since the electronic wavefunction of the edge states is 
more and more confined to an edge of a ZGNR when 
its wavevector k approaches the zone boundary, it pro- 



FIG. 3: (Color online) Calculated band structure and energy 
gap of ZGNRs. (a) The LSDA band structure of a 12-ZGNR. 
The up and down spin states are degenerated for all the bands, 
and the top of the valence band is set at zero. The symbols, 
A and A 1 denote the direct band gap and the energy gap at 
the zone boundary, (b) Variation of direct band gap with the 
width of ZGNRs. The open symbols denote the LSDA results 
while the solid symbols are the GW results, (c) Variation 
of the energy gap at the zone boundary with the width of 
ZGNRs. The symbols have the same meaning as those in (b). 



vides a possibility to see how the self-energy correction 
evolves with the localization of the electronic state. We 
plot the charge distributions of three electronic states of 
the first conduction band with different wavevector k and 
their corresponding self-energy correction values defined 
as E QP -E LSDA in Fig.gl It is clear that the self-energy 
correction is enhanced when the state is confined at the 
edge as shown from Fig. 0] (b) to (d). Because of the 
1/r nature of the Coulomb interaction, the self-energy of 
a state is sensitive to the localization of the wavefunc- 
tion. Therefore a larger self-energy correction is found 
to the more localized edge state. As a consequence, the 
dependence of the GW correction on wavevector signifi- 
cantly changes the band dispersion in ZGNRs from that 
of LSDA calculations. A smaller effective mass and bet- 
ter mobility for the carriers are expected in ZGNRs for 
the GW bands as compared to the LSDA ones. 

Recently, several experiments related to the quasiparti- 
cle band gap in GNRs have been reported 22j, [23| • They 
have not only proven the existence of finite band gap in 
GNRs but also shown a larger gap when the width of the 
GNR decreases. Within a range of width of GNRs of 15 
nm to 90 nm, a E g oc 1 /w relation is observed. This find- 
ing agrees qualitatively with our GW results. However, 
the experimental data are for the wider GNRs where the 
widths are far from the range of widths of our calcu- 
lated GNRs (0.4-2.4 nm). In addition, all the GNRs in 
the experimental case are etched by the oxygen plasma, 
which could be different from our hydrogen passivated 
GNRs. Therefore it is difficult to compare our GW re- 
sults with current experimental data directly. On the 
other hand, considering that the origin of the enhance- 
ment of the self-energy correction in GNRs is the quasi- 
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FIG. 4: (Color online) (a) Variation of GW correction (the 
difference between the quasiparticle gap and the LSDA gap) 
with wavevector of electronic states (k = 0.375, 0.4375 and 
0.50 in units of 2tt/c) in an 8-ZGNR. (b), (c) and (d) are 
the charge distributions of the conduction state with the cor- 
responding wavevectors in (a). Because the up and down 
spin components are degenerated but localized on a different 
sublattice, we plot the charge distribution of only one spin 
component. 



one-dimensional geometry and weakened screening, we 
expect that other passivating atoms or molecules do not 
change the physics here significantly. With advance in 
experimental techniques, it is very possible that smaller- 
sized and hydrogen-passivated GNRs will soon be fabri- 
cated. A comparison between our first-principles results 
and experimental data can then be made. 

In conclusion, we have performed a first-principles 
Green's function calculation within the GW approxima- 
tion to obtain the quasiparticle band gaps in GNRs. Due 
to the enhanced electron-electron interaction in these 
quasi-one-dimensional systems, a significant self-energy 
correction is found for both armchair and zigzag GNRs. 
The quasiparticle energy of states near the band gap in 
ZGNRs is found to be wavevector sensitive, and this gives 
rise to a larger band width and smaller effective mass for 
carriers in ZGNRs. The calculated quasiparticle band 
gaps are within the most interesting range (1-3 eV for 2- 
1 nm GNRs) and give promise for applications of GNRs 
in nanoelectronics. 
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